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Abstract 

Modern Computational Fluid Dynamics (CFD) 
techniques were used to compute aerodynamic forces 
and moments of the Space Shuttle Orbiter in specific 
portions of contingency abort trajectory space. The 
trajectory space covers a Mach number range of 3.5-15, 
an angle-of-attack range of 20°-60°, an altitude range 
of 100-190 kft, and several different settings of the 
control surfaces (elevons, body flap, and speed brake). 
Presented here are details of the methodology and 
comparisons of computed aerodynamic coefficients 
against the values in the current Orbiter Operational 
Aerodynamic Data Book (OADB). While 
approximately 40 cases have been computed, only a 
sampling of the results is provided here. The computed 
results, in general, are in good agreement with the 
OADB data (i.e., within the uncertainty bands) for 
almost all the cases. However, in a limited number of 
high angle-of-attack cases (at Mach 15), there are 
significant differences between the computed results, 
especially the vehicle pitching moment, and the OADB 
data. A preliminary analysis of the data from the CFD 
simulations at Mach 15 shows that these differences can 
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be attributed to real-gas/Mach number effects. The 
aerodynamic coefficients and detailed surface pressure 
distributions of the present simulations are being used 
by the Shuttle Program in the evaluation of the 
capabilities of the Orbiter in contingency abort 
scenarios. 


Nomenclature 

C p = Pressure coefficient 

C A = Axial force coefficient 

C HBF = Hinge moment coefficient (body flap) 

Chei = Hinge moment coefficient (inboard eleven) 
C HE0 = Hinge moment coefficient (outboard elevon) 
C Hsn = Hinge moment coefficient (speed brake) 

C m = Pitching moment coefficient 

C N = Normal force coefficient 

L b = Body length (ft) 

M = Mach number 

Pr = Prandtl number 

q = Dynamic pressure (psf) 

Re = Unit Reynolds number (1 /ft) 

Sc ~ Schmidt number 

7* = Static temperature (°R) 

V = Velocity (ft/sec) 

Z = Altitude (ft) 

x = Axial coordinate (ft) 

y = Span coordinate (ft) 

z = Normal coordinate (ft) 

a = Angle of attack (°) 

Y = Ratio of specific heats 

AC m = Pitching moment coefficient increment 

S BF = Body flap deflection (°) 

<5 f = Elevon deflection (°) 

&sb = Speed brake deflection (°) 

<p - Meridional angle (°) 

p = Mass density (slug/ft 3 ) 

Subscripts 

t = Turbulent 

sc = Freestream 
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Intro ductio n 

Assessment of the Space Shuttle Orbiter 
performance during potential contingency abort 
scenarios is an essential element of future upgrades to 
improve the safety and reliability of the Shuttle fleet. 
Improving the Shuttle contingency performance is an 
important part of the overall Space Shuttle Orbiter 
safety improvement plan. 1 Orbiter abort scenarios are 
currently being re-examined to improve safety for 
Return to Launch Site (RTLS) aborts (see Fig. 1). The 
current two- and three-engine out trajectories for RTLS 
aborts include regions (black zones) where loss of the 
Orbiter and crew is predicted due to vehicle limitations 
and constraints (see Fig. 2). The opportunity to recover 
the crew and vehicle in these abort situations could be 
improved if the operational capabilities of the Orbiter 
outside of its nominal range were better understood. 
Therefore, a key to minimizing the black zones is to 
enhance our understanding of the aerodynamic 
capabilities of the Orbiter by providing better estimates 
of the aerodynamic loads the Orbiter would experience 
during abort scenarios. 

Many design assumptions that were intended to 
provide the required margin of safety prior to the first 
Shuttle flights have been reassessed with improved 
tools and flight data. In addition, the performance limits 
of the Orbiter are now well established for nominal 
operating conditions. The aerodynamic environment of 
the current Orbiter is documented in the Orbiter 
Aerodynamic Data Book (OADB). The OADB is 
primarily based on wind tunnel measurements 
extrapolated to flight conditions and updated with flight 
data where available. Much work was undertaken to 
validate the design aerodynamic database (i.e., the pre- 
flight aerodynamic database) with measurements from 
flights STS-1 through STS-5. After the first flight, the 
real-gas influence on the aerodynamic characteristics of 
the Orbiter at high Mach numbers was investigated 
fully including computational efforts of several 
researchers* that accurately captured the flow physics in 
the nominal flight regime. The validity of the OADB is 
not in question for nominal flight conditions. The 
results from the present work, however, suggest that 
there are selected areas where the current OADB could 
be improved in off-nominal conditions, where the 
current data-book values have been extrapolated from 
wind tunnel and flight data. CFD has matured 
substantially since the original database was 
constructed and can now be used to evaluate and update 
the existing aerodynamic loads database. 

Early RTLS trajectory assessment for contingency 
abort improvement requires aerodynamic loads in the 
freestream Mach number (AfJ range of 6.0 to 15.0 and 


angles of attack (a) ranging from 20° to 60°. The 
current loads (surface pressure distributions) database 
does not cover all of these conditions and is mostly 
limited to Mach 3.5 and below. 

In order to establish the expanded aerodynamic 
database, full Navier-Stokes flow solvers were 
employed to predict aerodynamic loads on the Orbiter 
and its control surfaces. The scope of this paper is to 
outline the process used to compute and verify the CFD 
solutions that cover the required Mach number and 
angle-of-attack ranges for abort scenarios with a variety 
of control surface deflections. In addition, we show 
representative results for some of the cases and discuss 
reasons for occasional differences between current 
OADB data and the computations. A total of 31 cases 
were identified to cover Mach numbers from 3.5 to 
15.0, four angles of attack (20° to 60°), three 
elevon/body flap combinations, and two speed-brake 
deflections. In addition, seven benchmark check cases 
were originally chosen to compare CFD results with 
available wind tunnel and flight data. An eighth 
benchmark case was added later. 

Computation of the production aerodynamic loads 
was split between Boeing Reusable Space Systems 
(BRSS) and NASA Ames Research Center (ARC) 
while NASA Johnson Space Center (JSC) coordinated 
the effort. The Mach 3.5, 10.0, and 15.0 conditions 
were primarily a responsibility of ARC. The overall 
strategy implemented by ARC to compute the 
contingency abort flight conditions uses methodologies 
derived from the X-33 program/ 1 " The ARC team 
developed a set of automation tools to facilitate setting 
up cases, convergence and quality checks, and surface 
pressure and hinge moment data comparisons. 

All results were presented to a review board 
consisting of experts in aerodynamics, structural 
analysis. Space Shuttle operations and CFD. This board 
(the Orbiter Aerodynamics Panel) reviewed the entire 
process to determine if the data was of sufficient quality 
to be released to the Orbiter structural analysts. 

Meth o d o log y 

The matrix of cases chosen to generate the 
necessary data, including the benchmark cases, is 
shown in Table 1. The benchmark cases include one 
wind tunnel case (OA-146) at Mach 3.5 and angle of 
attack of 18.5°, three flight cases at Mach 3.5, 6.0, and 
15.0, and three additional flight cases at Mach 15.0 to 
complete an angle-of-attack sweep from 30° to 60°. 
Although integrated aerodynamic forces and moments 
are available for the flight cases, only the comparison 
against the wind tunnel case will be presented in this 
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paper. This case is critical because it also allows for 
surface pressure comparisons between the CFD and the 
experimentally measured values at data port locations 
on the wind tunnel model. 

The production run matrix spans critical regimes of 
the contingency abort trajectories in the range of Mach 
numbers from 3.5 to 15.0 and angles of attack from 20° 
to 60°. The Mach 3.5 and 6.0 cases include a speed 
brake deflected at 69.76°. The speed brake is not 
deflected for the higher Mach number cases. The matrix 
also includes negative and positive control surface 
deflections for both the elevons and the body flap. The 
Mach 3.5, 10.0, and 15.0 runs and a limited set of Mach 
6.0 runs were calculated using the General 
Aerodynamic Simulation Program ( GASP ’). 

In order to perform the simulations efficiently 
within a reasonable time and to assure the quality of the 
computations, an integrated process was developed 
with numerous check points to identify, capture and 
correct any potential mistakes. This process includes 
grid generation, CFD simulations and also the 
translation of the CFD data sets into data formats that 
can be readily used by the structural evaluation group. 
The overall approach is outlined in Fig. 3 and consists 
of six phases: geometry certification and grid 
generation; CFD computations; post processing and 
analysis; conversion from CFD to Model Modal Loads 
and Structures (MMLS) format; aero analysis of 
MMLS formatted data sets; and the Orbiter 
Aerodynamics Panel review and approval. It is a 
requirement that the surface grid be certified as 
conforming to the Orbiter CAD definition provided by 
JSC. The initial grid generation phase includes the 
construction of surface and volume grids that 
topologically model the geometric details of the Shuttle 
Orbiter, including the control surfaces. A number of 
tools are available for analysis and verification. All 
computations are quality checked and cleared for post- 
processing. Tools were created as necessary to bridge 
the gap between the CFD results and the formats 
required for structural analysis. The converted results 
are then used for aero analysis before review, approval 
and release to the structures group. Details of the 
general approach are discussed in the following 
sections. 

Geometry Verification and Simplification 

Grid generation is a time consuming process. The 
overall strategy implemented here addresses all the grid 
requirements outlined in the production run matrix 
(Table 1) for combinations of control surface 
deflections within the schedule constraints. The true 
surface representation for the Orbiter was provided by 


JSC in CAD/IGES format. Individual surfaces were 
then stitched together and trimmed using Hyper Mesh 1 ' 1 
to provide a single “water tight” surface definition for 
the grid generation process. The resulting surface was 
then triangulated, capturing all of the geometric details, 
including control surfaces, engines and nozzles in the 
base region. The requirements for the surface grid are - 
(1) conformity of the triangulation to the CAD surface, 
and (2) grid smoothness. The surface grid was 
compared against the IGES by JSC and the wing root 
section near the fuselage-to-wing glove-fairing surface 
was found to deviate from the true surface. This 
tolerance error was introduced during the trimming 
process of the wing surfaces against the fuselage 
surfaces. A corrective action was taken to locally refine 
the surface grid and project it on the CAD geometry. 
The surface grid was then judged to conform to the 
CAD geometry and was certified for CFD simulations 
by JSC. 

At this point, certain geometric simplifications 
were made that facilitated the grid generation process 
without sacrificing solution fidelity. Regions where the 
Orbiter geometry was approximated are highlighted in 
Fig. 4a, and detailed views are shown in Figs. 4h-e. The 
true geometry at the wing tip (Fig 4b) is characterized 
by a backwards-facing step. To resolve the local flow 
phenomena, a local topological substructure would 
have been required. The contribution of this local flow 
structure on the global aerodynamics was assumed to be 
negligible. An edge-fairing surface was used to 
suppress this geometric detail and locally idealize the 
geometry as shown in Fig. 4c. For the low Mach 
number cases (3.5 and 6.0) the speed brake is deflected 
which would expose a gap on both sides of the control 
surface. As shown in Fig. 4J a gap fairing filler surface 
is used to cover the gap regions, which permits the 
same topology to be used for both the undeflected and 
deflected speed brake. The surface grid accurately 
captures the speed brake thickness and curvatures on 
either side at the tip of the control surface. The last 
geometric approximation is highlighted in Fig. 4c. All 
the geometric details and the engine nozzles in the base 
region were simplified with a shroud cover. It is 
important to note that the windward, side, and trailing 
edge surfaces of the body flap are exactly preserved. In 
the base region the local pressure is expected to be 
orders of magnitude lower in comparison to the 
pressure on the windward side of the vehicle, which is 
the main contributor to aerodynamic forces and 
moments. Unpublished solutions obtained using the 
USA code with all the geometric details in the base 
region, are used to assess this geometric simplification. 
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Volume Grid Generation 

For volume grid generation, a strategy was adopted 
to provide grids for a large number of high-fidelity 
solutions that have differing grid requirements for 
control surface deflections and angles of attack. The 
grid generation tools used were GridPro n and 
Gridgen. 12 GridPro was used to build the complex 
topological structure and Gridgen was implemented for 
local grid modifications. The volume grid consists of 
three topological layers: the inner, intermediate, and 
outer. The inner layer shown in fig. 5a was designed to 
capture the boundary layer and contain the topological 
structure around all the control surfaces. Localized 
topologies were built in this layer to allow for 
deflection of the control surfaces by local substitution 
without manipulating the overall volume grid. In the 
inner layer, the normal clustering and the number of 
points can be adjusted as needed to capture the viscous 
gradients without altering the overall volume grid. The 
intermediate layer shown in Fig. 5b reflects the 
topological complexity due the surface geometry. Most 
of the complexity is in the aft part of the vehicle 
associated with the deflected speed brake and elevons. 
A strategy was developed to degenerate and redirect 
topological complexity away from the nose of the 
vehicle where regions of high grid density associated 
with the control surfaces should not propagate. The 
outer layer shown in Fig. 5c was designed to efficiently 
contain the bow shock for each angle-of-attack. 
Modifications to the outer layer grid, related to changes 
in angle of attack, do not disturb the other layers. Figure 
6;i highlights the surface grid smoothness on the wing, 
as well as the strategy of redirecting grid points from 
the inner elevon gap region to the outer elevon gap. 
Figure 6b shows the inner grid layer where the grid 
density normal to the surface and distributions could be 
adjusted accordingly to capture the viscous layer. 
Figure 6c highlights the localized topology where local 
volume grids for deflected control surfaces could be 
replaced. The final grid is comprised of 365 elementary 
blocks (see Figs. 7a, b) that are merged down to 38 
point-wise matching zones, and contains about 2.2 
million grid points. It was later determined that a wake 
grid was required for the high angle-of-attack (50° and 
60°) cases in order to contain the subsonic portion of 
the wake flow. This grid was built using Gridgen and 
added 1.4 million points to the volume grid. 

CFD Simulations 

The GASP code solves the full Navier-Stokes 
equations using finite volume spatial discretization. 
GASP includes a wide selection of physical models and 
numerical algorithms. Physical model selections are 
based on the free-stream conditions, and the choice of 


numerical algorithm is based on grid topology and 
sensitivity studies performed previously. 1 11 The ideal- 
gas model (y=1.4) is used for cases up to Mach 6.0. For 
higher Mach numbers, real-gas effects are important 
and are modeled using Park’s 1 five-species air model 
and finite-rate chemistry. The inviscid fluxes are 
computed using third-order Van Leer flux vector 
splitting with the min-mod limiter. The discretized 
system of equations is advanced to a steady state using 
two-factor approximate factorization. For all turbulent 
cases, the Baldwin-Lomax turbulence model 4 (with 
compressibility corrections') is used. The implemented 
physical models and numerical algorithms are 
organized in Tabic 2. GASP has been applied 
successfully in the past to predict aerothermal 
environments for the Shuttle Orbiter, X-33, X-34, X-38, 
Mars Pathfinder, Mars’03 and other planetary mission 
studies GASP has been a very useful tool in relating 
facility measurements to flight and in providing insights 
into viscous, real-gas, and Mach number effects. In 
addition, GASP simulations performed for wind tunnel 
conditions (A/*=6.0 and 20.0) for the X-38 vehicle 
predicted aerodynamics well. 5 

Post-Processing 

Global aerodynamics were extracted from the 
computed solutions by integrating pressure and shear 
stress over the Orbiter surface. Certain quantities, 
notably pitching moment in a near-trim attitude, are 
extremely sensitive to interpolation errors and other 
biases during extraction of the aerodynamic forces and 
moments, and care must be taken to eliminate these 
errors from the integration. For example, early in the 
project there was a one-inch error (less than 0.1% of 
vehicle length) in the stream-wise location used for the 
moment reference center. This changed the computed 
pitching moment by as much as 12% near trim 
conditions. Control surface aerodynamics were 
extracted in a similar manner, with the integration 
occurring over only that subset of the surface grid that 
defines the control surface. 

Historically, a stumbling block in the 
communication between the CFD and structural design 
communities has been the difficulty of transforming 
CFD solution data into a format appropriate for 
practical vehicle design. In this project, a format that 
could be used by the MMLS suite of codes was 
requested by the Aero and Structures groups. Values of 
surface pressure coefficient ( C p ) were to be provided at 
the fixed locations (data ports) on the Orbiter surface. 
For design requirements the vehicle data ports are 
grouped into four segments: the fuselage; upper and 
lower wings; and the vertical tail; with each referenced 
to their own local coordinate system. The provided 
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locations of the data ports were mapped onto the actual 
Orbiter surface, and the resulting locations were 
verified by the structures team. Although there were 
some initial problems with locating points appropriately 
on control surface leading and trailing edges, these 
issues were resolved after a few iterations between the 
CFD and structures teams. The CFD data were then 
interpolated to the final data-port locations using an in- 
house tool, and the resulting data were compared to the 
original CFD solutions to ensure that the interpolation 
was accurate. 

R ggu)tg 

Qrid Re so lu tion 

All CFD simulations using GASP were performed 
on two levels of grid refinement. The grid sequencing 
strategy was limited by the multi-block structure. The 
complex topology necessitated coarsening the grid in 
all three coordinate directions simultaneously. As a 
result, the coarse and fine grids differ by a factor of 
eight in grid density. The coarse grid level (Fig. 8a) 
contains 275 thousand points and the fine grid (Fig. 8h) 
has 2.2 million points. Figure 8c shows computed 
surface pressure coefficient on both grids for the Mach 
10, a=40° case. From the figure we see that both grids 
produce very similar C p distributions. Also shown in 
Fig. He are the integrated axial (C*) and normal (C^) 
force and pitching moment ( C m ) coefficients on the two 
sequence levels. Differences in the normal force and 
pitching moment coefficients are very small. The 
largest difference is observed in the computed axial 
force. This result is not surprising since the coarse grid 
has half the resolution in the surface-normal direction. 
As a result, the viscous contributions to the axial force 
are over predicted on the coarse grid. 

As a standard practice, all the CFD solutions are 
checked for steady-state convergence. In addition to 
continuously monitoring the residual, the integrated 
loads and moments are checked at regular intervals. 
The solution is converged when the integrated loads 
and moments reach asymptote values and no longer 
change. A more detailed convergence check is 
performed by monitoring changes in surface C p 
distribution. Differences are taken between calculated 
surface C p distributions every four hundred iterations 
and images are generated to look for asymptotic 
convergence. This process is illustrated in Fig. 9 for the 
Mach 3.5, a- 18.5° benchmark case. The maximum 
difference in C p on the wind and lee surfaces between 
solutions four hundred iterations apart is less than 
0 . 001 . 


Base Contributions 

The impact of the base geometry approximation on 
aerodynamic forces and moments is addressed here. 
Figure 1 0a shows a schematic of the ARC 

approximated base geometry. All the geometric details 
in the base region, including the engine nozzles, were 
replaced by a shroud cover. The justification for this 
simplification is that the base region is exposed to 
pressure levels that are orders of magnitude lower than 
the windward pressures. Therefore, the wind-side 
pressure distribution is the main contributor to the 
integrated loads and moment coefficients. In order to 
quantify the effect of this simplification, base 
corrections to all global aerodynamic quantities were 
extracted from USA solutions (computed by the BRSS 
team; not discussed in this paper), which included all of 
the geometric details in the base region. The effect of 
the base simplification is discussed here only in terms 
of the integrated moment coefficient. The effect on 
other aerodynamic coefficients is much smaller. Figures 
Hike show pitching moment as a function of angle of 
attack for a constant Mach number including the shuttle 
data-book values (solid lines) along with their reported 
uncertainty (dashed lines). The solid circles represent 
the pitching moment values that exclude the 
contribution from the approximate base geometry, i.e., 
the surface aft of the solid line in Fig. I()a. The solid 
triangles represent the pitching moment values that 
include the contribution from the approximate base 
geometry. The open circles are the GASP results that 
have been “corrected” to include a base contribution 
taken from the corresponding USA solution. In Fig. 

I Oh, comparisons are shown for the Mach 3.5 case, 
since the base contribution is anticipated to be larger at 
lower Mach numbers. Figure 10b shows that the 
differences are small and that the base correction tends 
to push the computed pitching moment toward the 
OADB values. Figure 10c shows the same comparison 
for the highest Mach number case. In this case, the base 
contribution to the overall integrated moment is 
negligible. This assessment re-affirms that the 
geometric simplifications in the base region do not 
sacrifice solution quality. It is important to clarify here 
that although the base geometry was simplified, the 
high angle-of-attack cases (50° and 60°) were computed 
with a wake grid. A wake grid was necessary for these 
cases to provide the appropriate boundary conditions to 
capture the flow physics as the flow expands rapidly 
from the windward side into the base region. For all the 
low angle-of-attack cases (20°-40°) a wake grid was not 
included in the computation. For these cases it was 
effectively assumed that the base contribution to 
aerodynamic forces and moments is zero. 
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Comparison of the individual cases with the OADB 
provided a basis to look at the CFD simulations more 
closely. For most cases, the solutions were in good 
agreement with the OADB, and were deemed 
acceptable once they had passed the internal quality 
review. In those cases where there were observable 
differences, the CFD simulations were First re-evaluated 
and then re-analyzed to understand the cause of the 
differences. In addition, the origin and makeup of the 
OADB for these cases were reviewed to determine the 
fidelity of the data book for nominal and off-nominal 
conditions. In the following sections we present sample 
results and focus on a flight regime where CFD predicts 
interesting results. 

Mach 3 ,5 Computation s 

All Mach 3.5 cases, except for the -15° (trailing 
edge up) elevon deflection cases, shown in the 
production run matrix (Table 1) were computed 
successfully using the GASP code and the results were 
compared against the data-book values. The Mach 3.5, 
d E = -15° cases were not computed due to schedule 
constraints. Instead an existing volume grid for the -10° 
elevon deflections was used to compute the Mach 3.5 
schedule B (40° — 60°) cases. An extrapolation scheme 
has been established based on experimental data and 
CFD (<5 £ = -10°) to provide integrated forces and 
moments as required for those cases. A representative 
set of cases for the Mach 3.5 schedule B angle-of-attack 
sweep is showrn in Fig. 14. For the cases presented, both 
the elevon and body flap control surfaces are 
undeflected. The speed brake is deployed to 69.76°. 
The comparisons are very good in general. Trends 
follow the data book and comparisons are on or within 
the uncertainty bands for the normal {C N ) and axial (C A ) 
forces, pitching moment (C m ) coefficients and control 
surface hinge moments, inner elevon (C HEl ), outer 
elevon ( C HE0 ), and body flap ( C HBF ). Tabic 4 
summarizes the comparisons of all force and moment 
coefficients including hinge moments for all control 
surfaces for the cases shown in Fig. 14. In Table 4 the 
moment coefficient for the Mach 3.5, a= 50° with 
undeflected elevons and body flap and the a = 50°, S E = 
-10° cases are marked with marginal acceptability but 
are sufficiently close to the data-book values for the 
contingency abort loads assessment purposes. The same 
classification was also given to the Mach 3.5, a = 60°, 
d E = -10°, undeflected body flap (& BF - 0°), elevon 
hinge moment comparisons. Significant differences are 
observed for the body flap hinge moments at the a = 
50° and 60°, S BF —11.7° cases. For these cases the 
global force and moment comparisons are very good 
and local differences in the body flap region are related 
to strong windward side sonic flow at high angles of 


attack. Detailed investigations are currently underway 
to understand the precise nature of these differences. 

Mach 6 Computations 

The forces and moments at Mach 6 are compared 
to the data-book values in Fig. 15. The cases presented 
here correspond to undeflected elevons and body flap. 
The speed brake is deflected to 69.76°. The angle-of- 
attack sweep shown corresponds to schedule C 
(20°-60°). The overall comparisons between the 
computations and the OADB values are excellent. The 
CFD computations capture the pitching moment vs. 
angle-of-attack slope change (cusp) that occurs at 50°, 
consistent with the OADB measured values. For these 
high angle-of-attack conditions the flow in large 
portions of the windward side shock layer becomes 
subsonic and as a result the pressure distribution 
changes. This significantly changes aerodynamic 
characteristics of the Orbiter, which is evident in the 
integrated aerodynamics by the cusp shown in Fig. 15c. 
More details on this phenomenon will be presented later 
in this paper. The comparison summary is shown in 
Tabic 5 All Mach 6.0 data were accepted by the 
Orbiter Aerodynamics Panel, including the speed brake 
hinge moment value at a = 20° where the data were 
identified as marginally acceptable for assessment 
purposes. 

Mach 10 Computations 

A single case at Mach 10.0 was included in the 
production run matrix to provide an anchor point at a = 
40° where flight data exist. This check case was 
computed using GASP to compare with OADB 
integrated forces and moments and to verify Mach 
trends. The results are shown in Figs. 16a -g. The 
integrated aerodynamics were found to be in good 
agreement with the data book. CFD results were 
reasonable both in magnitude and Mach trends. 

15 CppipptatiQns 

In order to assess the sensitivity of the computed 
solutions to the choice of transport models at Mach 15, 
turbulent flow simulations were performed at a = 30° 
and 40° in addition to laminar simulation for all angles 
of attack. The Baldwin-Lomax turbulence model ' with 
compressibility corrections 4 was used in the 
computations. Since GASP does not currently possess 
the ability to model transition, the flow was assumed to 
be fully turbulent from the tip of the nose. Tabic 6 
shows the computed aerodynamic coefficients for the 
two angles of attack. Compared to the corresponding 
aerodynamic coefficients from the laminar 
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computations, turbulent predictions show substantial 
increases in the axial force coefficients and small 
changes in the normal force coefficients. Further, the 
small changes in pitching moments are “stabilizing” 
(i.e., increasingly negative compared to the OADB 
mean values). The observed changes are not wholly 
unexpected. The assumption of fully developed 
turbulence at the nose of the vehicle is questionable. It 
must be borne in mind that at Mach 15, the unit 
Reynolds number (based on free stream quantities) is 
small (-34,000/ft) and turbulence is not expected to 
develop as far forward as the nose of the vehicle. 
Transition would occur, if at all, towards the aft of the 
vehicle. Given that the center of gravity of the vehicle 
is located far aft, the moment contributions from the 
nose region of the vehicle can be very large due to the 
large moment arm. Any significant increases in the 
shear load in this region will lead to large increases in 
the computed moment. Any increases aft would 
probably have a smaller impact due to the shorter 
moment arm to the center of gravity of the vehicle. The 
assumption of a laminar flow is bolstered by the fact 
that at a = 40°, the predicted force and moment 
coefficients are in good agreement with the OADB 
values which are extracted from the Shuttle flights. 

In addition to the turbulent results, Tabic 6 shows 
the contributions of the non-viscous (pressure only) 
terms to the aerodynamic forces and moments for the 
laminar cases. It is evident from the tabulated results 
that the viscous contributions to the forces and 
moments are very small with the sole exception of the 
axial force. The contribution of the wall-bounded shear 
layer to the axial force is significant as expected. 

Based on these sensitivity studies, all of the Mach 
15 cases were computed assuming the flow to be 
laminar. The additional discussion of the Mach 15 cases 
that follows is based on the computed laminar solutions 
(either ideal-gas or real-gas model). 

CFD to OADB comparisons of results of select 
cases from the Mach 15 production run matrix (shown 
in Table 1) are presented in Figs. 17a-f. The integrated 
loads and moments are plotted against the schedule A 
(30° — 60°) angle-of-attack range. The comparisons are 
very good in general, and the computed values are on or 
within the uncertainty bands. The computed data follow 
the OADB trends except for the integrated pitching 
moment coefficient (C m ) at a=60°. For this particular 
case, GASP predicts a more negative (nose down) 
pitching moment and is divergent from the OADB 
value. The OADB trend shows the cusp that was seen 
earlier in the Mach 3.5 and Mach 6 data. For the Mach 
3.5 and 6 computations, the cusp was captured by the 
CFD and comparisons with the OADB data were very 


good. Therefore, it was puzzling that the CFD seemed 
unable to capture the similar trend at the higher Mach 
number. Further studies were undertaken to resolve 
these differences between the computed results and the 
OADB values and the results of these studies are 
presented and discussed in the next section. 

The comparison summary for all the Mach 15 
cases is shown in Table 7. All Mach 15 data were 
accepted by the Orbiter Aerodynamics Panel, including 
the a= 60° pitching moment values discussed below. 
High Angle of Attack Comparison at Mach 15 

In the previous section, the results of the Mach 15 
computations were presented (see Figs. 17a-f)- The 
computed results, with the exceptions of the vehicle 
normal force and pitching moment coefficients at an 
angle of attack of 60°, compare favorably (i.e., within 
the uncertainty bands) with the OADB. At a = 60°, 
however, the computed pitching moment is found to be 
more “stabilizing” (i.e., more negative or nose down) as 
compared to the OADB. Further, the cusp (or a pitching 
moment reversal) in the OADB pitching moment curve 
(see Fig. i 7c) beyond a - 50° is not reproduced by the 
computed results. The large negative value of the 
predicted pitching moment at a - 60° means that a 
large restorative (nose up) aerodynamic force would be 
required in order to trim the vehicle at this angle of 
attack. This restorative force is much larger than that 
required using the OADB value of the pitching 
moment. It is interesting to note that the results of 
computations performed at Mach 3.5 and 6 agree very 
well with the OADB values and, more importantly, 
show the same trends as the OADB up to and including 
a = 60°. 

In Figs 18a-c, the pitching moment coefficient is 
plotted as a function of Mach number for angles of 
attack of 40°, 50° and 60°, respectively. At a - 40° 
(Fig. 18a), the pitching moment trend is in excellent 
agreement with that of the OADB. Again at a= 50° 
(Fig. 18b) the trend is in good agreement with the 
OADB. At a = 60° (Fig 18c), however, the computed 
pitching moment trend in Mach number is opposite to 
that of the OADB beyond Mach 6. The results of a 
preliminary investigation into resolving this 
aerodynamic behavior at a = 60° is now presented. 

We begin by examining the OADB data at high 
Mach numbers and angles of attack. The OADB mainly 
consists of pre-flight aerodynamic data gathered from 
extensive wind tunnel testing from either ideal-gas (air) 
facilities (for Mach 3.5 and 6.0) or inert-gas (He) 
facilities (for Mach 18-20). These basic data are 
supplemented with - (a) results of best available theory 
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(to account for, among others, viscous, Mach number, 
high-altitude, and real-gas effects), and (b) data from 
the various Orbiter flights. The supplemental data are 
added as increments to the basic pre-flight (ideal-gas) 
values. For example, the pitching moment coefficient 
can be written as a basic term (C m . flA5/c ) plus increments 
due to various effects, i.e., 

^rn-OADB ~ BASIC ^ AC ' m -Reul Gas A C m _ Viscous fnte rations 

A Altitude + '•* 

Cm-BAsic * s extracted from ground-based measurements 
and extrapolated across the Mach number range. This 
extrapolation leads to large uncertainties in the OADB 
data. Note that any narrowing of the uncertainty bands 
in the aerodynamic coefficients reflects the use of flight 
data in the OADB. The pitching moment increment due 
to real-gas effects (AC m _ Re(ll Cas as used in the OADB) at 
Mach 15 for a= 40° is +0.022 and is in fairly good 
agreement with that predicted using the GASP code 
(+0.0248) at the same angle of attack. 

It is important to note that - (1) no flight data exist 
at 50° and 60° angles of attack (however such angles do 
occur at points on the contingency abort candidate 
trajectories), and (2) the Orbiter is currently certified to 
fly at angles of attack less than 44° at Mach 15. 

The reversal in pitching moment at high angles of 
attack occurs when the shock-layer becomes subsonic 
over a large portion of the windward side of the vehicle. 
In Fig. 10 the sonic lines for the Mach 3.5 cases (for 
a= 40°, 50°, and 60°) are shown for the pitch plane of 
the vehicle. The sonic line attachment point moves aft 
with increasing angle of attack and at the highest angle 
of attack (60°), the sonic line anchors close to the 
trailing edge of the body flap. Further, at this angle of 
attack, the flow on the windward side of the vehicle is 
almost entirely subsonic as is expected. The CFD 
simulations at this Mach number match the behavior of 
the pitching moment data in the OADB. A similar pitch 
plane plot is shown for Mach 6 in Fig. 20. Again, the 
rearward movement of the sonic line is seen with 
increasing angle of attack. Further, the CFD simulations 
duplicate the pitching moment reversal at this Mach 
number. Note that all of the Mach 3.5 and 6 cases were 
computed using an ideal-gas model. Recall that the 
Mach 15 cases were computed using a real-gas model. 
To better understand the nature of real-gas effects at 
high angles of attack, the ideal-gas model was therefore 
used to compute additional solutions at Mach 15 for an 
angle-of-attack range of 30°-70°. In Fig. 21 the sonic 
lines for these Mach 15 cases are shown for the pitch 
plane of the vehicle. Again, the rearward movement of 
the sonic line attachment point is clearly seen. Further, 
at a= 60° and 70°, a substantial portion of the 


windward shock layer is subsonic. For a= 60°, the 
sonic lines for Mach 3.5, 6, and 15 ideal-gas cases are 
shown in Fig. 22. The impact of Mach number is 
clearly seen - the sonic line attachment point moves 
forward with increasing Mach number. In Fig. 23, the 
sonic lines in the pitch plane are shown for the real-gas 
cases at Mach 15. At a= 60° the sonic line does not 
attach as far aft on the vehicle as is observed in the 
ideal-gas computations. In the real-gas computations, 
the windward side of the vehicle develops a smaller 
subsonic zone (the sonic lines attach closer to the nose) 
and consequently, a pitching moment reversal is not 
observed in the computed results. The nature of the 
shock layer in the real-gas case is quite different from 
that in the ideal-gas case - the bow shock is closer to 
the body (thinner shock layer). The effect of the 
physical model (real vs. ideal) on the integrated 
pitching moment variation with angle of attack at Mach 
15 is shown in Fig 24. 

Mac h 15 Con cl q , s igns 

The real-gas predictions of the pitching moment 
are within the uncertainty band of the OADB up to an 
angle of attack of 50° with excellent agreement with 
flight data at a = 40°. While the real-gas computations 
do not exhibit the pitching moment reversal at high 
angles of attack, results of computations based on the 
ideal-gas model do. A more detailed investigation of 
real-gas effects is currently underway. While the 
observations made above are limited to Mach 15, what 
is not clear is whether the same trends will hold across 
the entire Mach number range where real-gas effects 
are important. In addition to Mach number and real-gas 
effects, the freestream density (or equivalently, the 
altitude) has an important role. The results from the 
present work suggest that the OADB could be improved 
or enhanced in the hypersonic, high angle-of-attack 
regime, an important regime from the point of view of 
contingency aborts. 


The general comparison between the computed 
integrated forces, moment coefficients and hinge 
moments against the OADB is summarized in Tables 4, 
5, and 7 for the Mach 3.5, 6.0, and 15.0 sets of cases, 
respectively. The overall agreement is very good. The 
majority of the computed data is typically inside or near 
the OADB uncertainty bands, and follow similar trends. 
In general, the computed moment coefficients are 
consistently lower than the OADB values, the normal 
force coefficients consistently higher, and the axial 
force coefficients typically agree very well. The Mach 
15, a = 60" pitching moment characteristics have been 
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discussed and differences between the CFD and the 
OADB trends were justified. All of the CFD -gene rated 
airloads have been released for structural analysis. 

The Mach 10 and 15 computations and 
comparisons to OADB data clearly demonstrate that 
real-gas computational methods are essential for high 
Mach number aerodynamic predictions. The production 
run cases computed by ARC resulted in reasonable 
agreements and close predictions consistent to the 
OADB values and trends. This general agreement 
points to increased confidence in the CFD process to 
predict aerodynamics in regimes where limited 
experimental and flight data exist. It also provides 
critical insights into physical phenomena such as Mach 
number, real-gas and viscous effects. This confidence 
level in the CFD data allows for a deeper understanding 
of the physics involved and tighter integration between 
wind tunnel, flight, and numerical data. 

This work has provided the post-processing tools 
for translating CFD results into formats useful by other 
disciplines in aerodynamic and structural analysis. 

The final recommendation of the working group 
was to release the CFD data generated for abort 
assessment and further aerodynamic and structural 
analysis. The detailed CFD process and computed 
results presented here were extensively reviewed by the 
Orbiter Aerodynamics Panel. The quality of the data 
was recognized by the review board and approved for 
release. 
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Table 1 Benchmark and production airloads run matrix 


Benchmark Cases 


M, 

a 


&BF 

&SB 

# of Cases 

ARC 

3.5, STS-93 Fit. 

18.5 

0° 

+11° 

69.76° 

1 

1 

3.5, OA-146 

18.5 

0° 

+11° 

69.76° 

I 

1 

5.5, OA-146 

25.2 

0° 

+ 77° 

69.76 ° 

7 

7 

6, STS-93 Fit. 

28.5 

0° 

+ 11° 

69.76° 

1 

1 

15, OADB 

A a 

0° 

+ 11° 

0° 

4 

2 

Total number of benchmark cases 




n+n 

5 [+11 

Production Cases 


a 

&E 

&BF 

^SB 

# of Cases 

ARC 

3.5 

B h 

0° 

0° 

69.76° 

3 

3 

3.5 

B 

0° 

-11.7° 

69.76° 

3 

3 

3.5 

B 

-15 o (-10°) 

0° 

69.76° 

3 

3 

6 

C c 

0° 

0° 

69.76° 

5 

5 

6 

c 

0° 

10° 

69.76° 

5 


6 

c 

10° 

0° 

69.76° 

5 


6 

c 

-10° 

0° 

69.76° 

5 


10 

40 

0° 

0° 

0° 

1 

1 

15 

A 

0° 

0° 

0° 

4 

4 

15 

A 

0° 

10° 

0° 

4 

4 

15 

A 

-10° 

0° 

0° 

4 

4 

15 

A 

10° 

0° 

0° 

4 

4 

Total number of production cases 




46 

31 

Total number of cases 




53f+U 

36/+ 7/ 


a A Schedule: a = 30°, 40°, 50°, and 60° 
h B Schedule: a = 20°, 30°, 40°, 50°, and 60’ 
c C Schedule: a = 40°, 50°, and 60° 


Table 2 Physical models and numerical algorithms implemented in GASP 



Ideal Gas Model 

5-species Model 

Thermodynamics 

II 

Eq. Stat. Mech. 

Laminar Transport 

Viscosity 

Sutherland’s Law 

Blottner’s Curve Fits 

Conductivity 

Pr - 0.72 

Eucken’s Relation 

Mass diffusion 

- 

Binary, Sc = 0.7 

Turbulent Transport 

Viscosity 

Baldwin-Lomax 

Baldwin-Lomax 

Conductivity 

(with compressibility) 
Pr = 0.9 

(with compressibility) 
Pr ( = 0.9 

Mass Diffusion 

- 

Sc, = 0.9 

Chemical Kinetics 

- 

Park Rates 

Inviscid Flux 

Van Leer 

Van Leer 

Accuracy 

3 "‘-Order (MUSCL) 

3 "‘-Order (MUSCL) 

Limiter 

Minmod 

Minmod 

Viscous Flux 

Central 

Central 

Accuracy 

2 nd -Order 

2 nJ -Order 

Wall BC 

Isothermal Wall 

Radiative Equilibrium 

Time Integration 

2-Factor AF 

(fully catalytic wall) 
2-Factor AF 
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Table 3 Input conditions for all cases in the production airloads matrix 



Z( ft) 

V* (ft/sec) 

q « (psf) 

A4 

p x (slugs/ft 3 ) 

T x ( n R) 

to (1/ft) 

Benchmark 

- 

2,223.9 

374.7 

3.5 

1.52E-04 

168.09 

2.5E+06 

Production 

100,000 

3,470.6 

199.8 

3.5 

3.32E-05 

409.36 

372,481 

Production 

135,000 

6,280.1 

129.5 

6 

6.56E-06 

456.11 

122,094 

Production 

140,000 

10,553.3 

295.3 

10 

5.28E-06 

463.68 

163,817 

Production 

190,000 

15,708.8 

95.88 

15 

7.77E-07 

454.73 

36,233 


Table 4 Mach 3.5 GASP to OADB summary comparison table 




Computed data within uncertainty bands of OADB 
Computed data slightly outside uncertainty bands of OADB 
Computed data well outside uncertainty bands of OADB 


Table 5 Mach 6 GASP to OADB summary comparison table 


&E < 5 BF &SB C s C m C A G HE j G heo C hbf G hsb 
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Table 6 Computed aerodynamic force and moment coefficients for Mach 15 


a 



C N 

C A 


Chej 

Cheo 

Chbf 

30° 

Pressure Only 

Laminar 

0.7519 

0.0499 


-0.0788 

-0.0454 

-0.1563 


Total 

Laminar 

0.7529 

0.0582 

0.0002 

-0.0790 

-0.0455 

-0.1567 


Total 

Turbulent 

0.7622 

0.0800 


-0.0870 

-0.0479 

-0.1591 

40° 

Pressure Only 

Laminar 

1.1757 

0.0541 

-0.0181 

-0.1483 

-0.0738 

-0.2681 


Total 

Laminar 

1.1770 

0.0630 


-0.1488 

-0.0741 

-0.2678 


Total 

Turbulent 

1.1874 

0.0885 

-0.0232 

-0.1555 

-0.0763 

-0.2720 


Table 7 Mach 15 GASP to OADB summary comparison table 



Computed data within uncertainty bands of OADB 
Computed data slightly outside uncertainty bands of OADB 
Computed data well outside uncertainty bands of OADB 
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Fig. 2 Current abort black zones for two- and three-engine failures during launch. 
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Wing Tip: Edge Fairing 


— * Base: Shrouded base 

Fig. 4 Geometric approximations. 
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Fig. 5 Multi-layer grid generation strategy. 



Fig. 6. Grid generation strategy: (a) surface grid smoothness and gap topology; 

(b) viscous inner-layer grid topology; and (c) replaceable control surfaces. 
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Fig. 7 ARC volume grid (2.2 million points without the wake grid - wake grid 

added 1.4 million points): (a) topological morphology, (b) 365 non-merged blocks. 



COARSE GRID 


(a) 


Wind Side 


CL - 0 .07*3, C, - 7. Cv - -0 




mm 


(c) r 


0.0455. C, • 1.2210. C„ . 0.0440 



0 OJ 0.4 O.a I U 1.4 I M U 7 


Fig. 8 Baseline grid fidelity verification: (a) wind-side comparisons of coarse (275K pts.) and 

fine (2.2M pts.) grid sequence levels; (b) lee-side comparisons; and (c) Mach 10 angle of 
attack 40° comparison of computed surface C p between coarse and fine grids. 
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Summary: 3450 Herat ione 

Max - 0.01 

Min - 0.00 


Wind I Lee 


Delta Cp between solutions 



- 0.030 - 0.025 - 0.019 -0014 - 0.008 - 0.003 0.003 0.008 0.014 0019 0.02 5 0 030 

Fig. 9 Example of steady-state surface pressure coefficient convergence monitoring 
(Mach 3.5 angle of attack 18.5° tunnel simulation.) 



Alpha (deg) Alpha (deg) 

Fig. 10 Impact of base geometry approximation on pitching moment: (a) schematic of the 
base shroud simplification; (b) BRSS base correction contributions to pitching 
moment at Mach 3.5; and (c) BRSS base correction contributions to pitching 
moment at Mach 15. 
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Fig. 11 GASP to wind tunnel experimental data wind centerline comparison at 
Mach 3.5 and angle of attack of 18.5°. 




(a) (b) 

Fig. 12 GASP to wind tunnel experimental data comparison at Mach 3.5 and angle of 
attack of 23.2°: (a) wind centerline comparison; and (b) circumferential cuts. 
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Fig. 14 Integrated load and moment comparisons between GASP and OADB at Mach 3.5, 
4=0°, $ BF =0°, 5^69.76°: (a) normal force coefficients; (b) axial force coefficients; 
(c) pitching moments; (d) inboard elevon hinge moments; (e) outboard elevon 
hinge moments; (0 body flap hinge moment; and (g) speed brake hinge moments. 
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Fig. 16. Integrated load and moment comparisons between GASP and OADB at 

Mach 10 (real gas), ^=0°, S 8F = 0°, ^=0°: (a) normal force coefficients; (b) axial 
force coefficients; (c) pitching moments; (d) inboard elevon hinge moments; 

(e) outboard elevon hinge moments; and (0 body flap hinge moment. 
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Fig. 17 Integrated load and moment comparisons between GASP, USA and OADB at 
Mach 15 (real gas), $£=0°, 5*^=0°, 5^0°: (a) normal force coefficients; 

(b) axial force coefficients; (c) pitching moments; (d) inboard elevon hinge 
moments; (e) outboard elevon hinge moments; and (0 body flap hinge moment. 
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